The Jamming point street-lamp in the world of granular media 
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The Jamming of soft spheres at zero temperature, the J-point, has been extensively studied both 
numerically and theoretically and can now be considered as a safe location in the space of models, 
where a street lamp has been lit up. However, a recent work by Ikeda et al [l[ reveals that, in the 
Temperature/Packing fraction parameter space, experiments on colloids are actually rather far away 
from the scaling regime illuminated by this lamp. Is it that the J-point has little to say about real 
system? What about granular media? Such a-thermal, frictional, systems are a-priori even further 
away from the idealized case of thermal soft spheres. 

In the past ten years, we have systematically investigated horizontally shaken grains in the vicinity 
of the Jamming transition. We discuss the above issue in the light of very recent experimental 
results. First, we demonstrate that the contact network exhibits a remarkable dynamics, with strong 
heterogeneities, which are maximum at a packing fraction ((>* , distinct and smaller than the packing 
fraction (l>\ where the average number of contact per particle starts to increase. The two cross- 
overs converge at point J in the zero mechanical excitation limit. Second, a careful analysis of the 
dynamics on time scales ranging from a minute fraction of the vibration cycle to several thousands 
of cycles allows us to map the behaviors of this shaken granular system onto those observed for 
thermal soft spheres and demonstrate that some light of the J-point street-lamp indeed reaches the 
granular universe. 

PACS numbers: 45.70.-n 83.80.Fg 



I. INTRODUCTION 



In a loose sense. Jamming describes everyday situa- 
tions where particles, objects, or people become dense, 
slow and rigid: one thinks of systems as different as 
sand piles, foams, or traffic jams as jammed systems 0. 
Significant progress was achieved in tlie field about ten 
years ago, when frictionless soft spheres at zero temper- 
ature were introduced as a minimal and seminal model 
for Jamming Q . This system has been extensively stud- 
ied and now serves as a point of reference for 
which Jamming has a precise meaning. Specifically, for 
models for which forces are represented by particle over- 
laps, the Jamming transition occurs when the system can 
only be compressed further by allowing overlaps between 
particles. From that point of view, it is essentially a 
matter of satisfying geometric constraints, and indeed, a 
formal identification with an algorithmic description has 
been established @, Q. For athermal systems, the Jam- 
ming transition is intrinsically out-of-equilibrium, and re- 
quires a precise characterization of the protocol used to 
prepare the system. However, many features of the tran- 
sition appear to be protocol independent [13], and for a 
given protocol on an infinite system, the Jamming tran- 
sition is entirely controlled by the packing fraction. The 
transition occurs at the so-called "point J", and coin- 
cides with the onset of isostaticity i.e., the number 
of steric and mechanical constraints imposed at the con- 
tacts exactly matches the number of degrees of freedom 
available to the particles. A number of geometrical and 



mechanical quantities exhibit clear scaling laws with the 
distance to point- J • One prominent signature of Jam- 
ming for systems of frictional particles is the singular 
behavior of the average number of contacts per particle 
z — zj (X {(j) — 4>j)°', where zj is equal to 2(i, where d is 
the space dimension, is the packing fraction at point 
J, and a ~ 0.5 0. The distribution of the gaps between 
particles displays a delta function at zero and a square 
root decay for increasing gaps, which is at the root of the 
singular behavior of the average contact number 0, [T2I - 

This framework has provided key physical insights into 
the nature of rigidity, and the structure/mechanics of 
disordered soft matter systems, such as emulsions [isj . 
foams [H,[I3 and grains 0, [l^. Of course, this ideal- 
ized model misses some of the key features of real sys- 
tems, such as friction for dry systems, interface effects for 
multiphase systems, or hydrodynamic interaction for sus- 
pensions. In particular, several works have shown that 
the Jamming scenario for static packings becomes more 
complex when friction comes into play |19l - [2lj . 

Furthermore, many systems of interest are not purely 
static: colloidal suspensions undergo thermal agitation; 
vibrated or flowing granular systems undergo mechanical 
agitation. Whether the Jamming framework is relevant 
in the presence of agitation remains an open, hotly de- 
bated issue d, H, [13] . On the one hand, one expects the 
singular nature of the Jamming point to be blurred. On 
the other hand, an anomalous dynamics is expected to 
occur, because particle motion, driven by agitation, may 
be influenced by the proximity of the singular point (see 
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FIG. 1: Temperature-Packing fraction phase diagram: 

At zero temperature, below Jamming, there is always a way 
to pack tlie particles without overlaps and the energy of the 
system is strictly zero. Above Jamming, there is no packing 
without overlaps and the energy, purely potential, is greater 
than zero. At finite temperature, the kinetic energy is never 
zero, and this feature blurs the picture. Contacts and overlaps 
are always present. 



figure [T]). A recent numerical study of harmonic splieres, 
in the presence of temperature, focuses on the dynamics 
in the region very close to the T = Jamming point [l| . 
The authors demonstrate that there is no singularity 
at finite temperature and identify a critical region in 
the vicinity of the Jamming point, where vibrational 
dynamics is maximally heterogeneous. They also re- 
port crossover lines, in the temperaturc-packing-fraction 
parameter space, between harmonic and non-harmonic 
regimes, originating at point J. Finally, on the basis of 
the dynamical behavior reported in the literature, they 
place existing colloidal experiments in the temperature- 
packing-fraction parameter space. Their main conclusion 
is that these experiments actually sit rather far from the 
critical regime of point J. 

In the past ten years, we have systematically inves- 
tigated horizontally shaken grains in the vicinity of the 
Jamming transition [23l - {26| . Starting with rigid brass 
disks, we observed very large heterogeneities of the dy- 
namics when focusing on minute dis plac ements on the 
order of 5 x 10^'^ grain diameters j23l - l25| : it was conjec- 
tured that these heterogeneities were connected to the 
dynamics at the contact scale. This was later confirmed 
using soft photo-elastic disks |26[. In the latter case, 
the signature of the dynamical heterogeneities was not 
as sharp, but we clearly demonstrated that the contact 
network exhibited a remarkable dynamics, with strong 
heterogeneities, which are maximum at a packing frac- 
tion (/)*, distinct and smaller than the packing fraction (t)\ 
where the average number of contact per particle started 
to increase. Furthermore, by varying the vibration fre- 
quency and observed that these two cross-overs merged 
in the zero mechanical excitation limit. 



The strong similarities shared by the above experimen- 



tal results with those reported in the numerical study 
of thermal soft spheres [l| call for further investigation. 
Indeed, one would like to know to the extent of over- 
lap between models of thermal harmonic spheres and the 
dynamical criticality of the granular packings. 

To address these questions, we present novel results 
spanning the short time (inner vibration cycle) dynamics 
of the photo-elastic soft disks using stroboscopic dynam- 
ics, and the longer times studied in previous studies [23l - 
[26j . In order to provide a background, we present a con- 
cise and reasonably complete picture of the dynamics, 
the forces and the contacts close to Jamming in the pres- 
ence of mechanical agitation. Within this context, we are 
able to: (i) conciliate hard and soft grain experiments, 
(ii) locate the granular experiment into a "temperature" - 
packing fraction phase diagram and, by so doing, discuss 
the relevance of the Jamming framework for describing 
granular systems. We conclude that our granular exper- 
iments do probe the same critical regime as those de- 
scribed by This, in turn, validates the use of soft 
sphere model to describe such systems close to Jamming. 

The paper is organized as follows. In section |TT1 we de- 
scribe the experimental set up in detail, emphasing the 
two modes of data acquisition, a fast one and a slower 
stroboscopic one, which allow us to explore the dynam- 
ics over six orders of magnitude in the timcscales. Sec- 
tion [ni] demonstrates that the force network is essentially 
isotropic and Section IIVI focuses on the dynamics of the 
contact network. This section summarizes the results 
already reported in [26| and supplements these results 
with the dynamical properties of the contacts at short 
timescales. Section|V]is devoted to the study of the mean 
square particle displacements. This study explicitly de- 
tails the data processing required to obtain a meaningful 
computation of these displacements. The quantitative 
results obtained in this section are the key elements of 
the discussion. Section IVII analyzes the dynamical het- 
erogeneities of the displacement field, relates them to 
those of the contact dynamics and show that they are 
embedded in the structural properties of the contact net- 
work. Finally, section IVIII synthesizes our observations, 
relates them to the previous study of brass disk exper- 
iments [23| - |25| performed in the same set-up, and dis- 
cusses the issue raised in the introduction, regarding the 
correspondence between thermal soft-sphere models and 
experiments on vibrated grains, in terms of dynamical 
behavior in the vicinity of point J. 



II. SETUP AND PROTOCOL 

We first review the details of the experimental set-up, 
which was adapted from [2^ in order to allow for the 
use of photo-elastic grains and the detection of contacts. 
We also review the different acquisition techniques, em- 
phasizing in particular, the fast image acquisition which 
allows us to characterize the dynamics within one vibra- 
tion cycle, as opposed to the previous studies, for which 
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FIG. 2: Sketch of experimental setup, (color online) 
(a): photo-elastic grains lighted by transmission by a polar- 
ized backlight, (b): confining cell, (c): wall piston, (d): 
force sensor, (e): micrometric stepper motor, (f): vibrating 
frame, (g): stepper motor ensuring vibration, (h): notched 
belt transmitting vibration, (i): shelf, (j): wall, (k): trans- 
lation stages. (1): stainless steel bars, (m): optical table, 
(n): CCD camera, (o): analyzers located on a rotating 
wheel, (p): shelf isolated from vibrations. 



one image per cycle was acquired in phase with the vi- 
bration. 



A. Setup 

The experimental setup is sketched in figure^ A bidis- 
perse mixture of ~ 8 000 4 mm and 5 mm photo-elastic 
disks (PSM-4) (a) lies on a glass sheet, and is confined 
in a cell (b), whose area can be tuned with a piston (c). 
The piston is attached to a force sensor (d) and a mi- 
crometric stepper motor (e). The packing fraction, 0, 
can be fine-tuned from 0.795 to 0.83, with a resolution 
oi S(j) = 5 X 10~^. Below the glass sheet, an LED back- 
light device, covered with a polarizing sheet, provides 
an intense, large, thin and uniform source of circularly 
polarized light. The glass sheet and the light are em- 
bedded in a frame (f), which vibrates horizontally with 
an amplitude a ~ 1 cm and frequencies / = 6.25, 7.5 
and 10 Hz. The oscillation is driven by a stepper mo- 
tor (g), a notched belt (h) and an eccentric revolving 
shaft, which are attached to a shelf (i), the stability of 
which is ensured by 300 kg of lead bricks ballast and a 



rigid bracket to the wall (j). The confining cell is me- 
chanically decoupled from the vibration devices. It is 
embedded in a larger frame, which in turn is attached 
to four manual micrometric translation stages (k). This 
ensures a precise leveling of the confining cell with re- 
spect to the oscillating board. The translation stages are 
attached to stainless steel bars (1), which are screwed to 
an optical table (m). Also attached to the optical table 
is a trigger. The trigger is made of a refiection photo- 
transistor/photo-diode device, together with a Schmitt 
trigger electronic circuit. The device is placed in front of 
the revolving shaft, where a piece of black tape has been 
taped; when the sensor is in front - respectively outside 
of - the tape, it delivers a 5 V, -respectively V signal. 
The phase of the trigger fall is chosen to be when the 
velocity of the plate is minimum and the Mark-to-Space 
ratio is adjusted in such a way that the transients of the 
stepper and the exposure times occur separately. 



B. Data acquisition 

We want to investigate the dynamics, both at short 
times, namely within the vibration cycles, and at long 
times, that is over several thousands cycles. Altogether, 
the experiment covers seven decades of time steps and, 
apart from the force sensor (d), all our data comes from 
image acquisition. We thus need to conduct two sepa- 
rate series of experiments, one with a fast camera, run- 
ning continuously, and one with a standard CCD camera, 
triggered by the motion of the oscillating plate. In both 
cases, we access both the position of the grains and the 
photo-elastic pattern inside the grains. This cannot be 
achieved simultaneously, and we need to adapt the ac- 
quisition schema in order to be as close as possible to 
this ideal situation. The camera is fixed on a shelf (p), 
lying on an optical table and isolated through a rubber 
gasket in order to reduce the transmission of vibrations 
and minimize blur on the pictures. 

To record the displacements and the force network dy- 
namics at short times, we use a fast camera (2000 frames 
per sec) with a resolution of 1024 x 1024 pixels, which 
record 1361 frames during up to 6 cycles of vibration, for 
the largest vibration frequency of 10 Hz. We successively 
acquire two movies, with and without introducing an an- 
alyzer in the field of view of the camera. Only a few tens 
of vibration cycles separate the two acquisitions. Since 
the dynamics is completely frozen (see below), the pack- 
ing barely moves, and synchronizing the two movies, we 
associate the photo-elastic pattern and the grains cap- 
tured on the white-light (no crossed polarizers) images. 

The long time dynamics is recorded with a high reso- 
lution (2048 X 2048) CCD camera (n) triggered in such 
a way that the images are taken in phase with the mo- 
tion of the oscillating board. Analyzers (o) located on 
a rotating wheel, with minimal inertia, are inserted in 
the field of view of the camera once every two cycles, 
using a triggered stepper motor, so that white-light (re- 
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FIG. 3: Interparticle force measurement, (color online), 
(a): Grain detection and photoelastic signal. The photoe- 
lastic signal is displayed using the following color code (blue: 
weak signal to red: strong signal). The grains are drawn in 
white, (b): Tesselation structure: Grains are drawn in black; 
Delaunay triangles are sketched in red; Voronoi vertices are 
linked together by blue lines; Voronoi vertices are linked to 
grain positions by green lines. Each contact is defined by four 
zones, a, b, c and d. (c): Contact force measurements: value 
of the measurement in each triangle. Redder colors corre- 
spond to higher forces and bluer colors to lower forces, (d): 
Contact forces: Redder colors and longer lines correspond to 
higher forces and bluer colors and shorter lines to lower forces. 



spectively cross-polarized) pictures are taken every odd 
(respectively even) vibration cycle. It is then straightfor- 
ward to match the photoelastic pattern to the white-light 
images of the grains. In order to minimize blur, the pic- 
tures are taken at the phase for which the board velocity 
is minimal, that is when it reverses direction. This is 
also when the acceleration is maximal. We shall see in 
the following that this has direct consequences on the 
contact number measurement. Also, the stepper motor 
that switches polarizers position is attached to the ceil- 
ing, to avoid the transmission of its vibrations to the 
camera. Finally, in order to prevent thermal expansion 
of the grains due to heating, the LED backlight is also 
triggered on the vibration and flashes only during 6 ms, 
which is also the time exposure of the camera. 

From the white-light images, we extract grain posi- 
tions, and diameters (black circles in figEl^b)), on which 
we perform Delaunay triangulation (red lines in figl3|b)) 
and Voronoi tesselation (blue lines in figlSKb)). The 
grain positions are obtained with a resolution of 0.5% 
of d. Once the grains have been detected, an estimate 
of the pressure within each grain is obtained by inte- 
grating the square gradient of the cross-polarized light 
intensity over the disc area. We denote by Gf this es- 
timate of the pressure in grain i. The resolution in 
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FIG. 4: Protocol and Calibration, (color online), (a): 

Calibration curve of the photoelastic signal, vs. piston 
force F/Mg. (b): Packing fraction, <j), and piston force, 
F/Mg vs. time during the overall experimental run (prepara- 
tion and acquisition). The vibration frequency is / = 10 Hz. 



each grain is not good enough to carry out a force in- 
verse algorithm for the photo-elastic problem and 
compute the forces at contacts. However, we can esti- 
mate them as follow. For each inter-particle contact, we 
use the two particle positions and the positions of their 
two common Voronoi vertices to build a patter of tri- 
angles, which we call a, 6, c and d (see figure [3l^b)). We 
then compute the spatial gradient of the associated cross- 
polarized image (see figure (Ha)), and we sum this signal 
within each of the triangles (see figure [3l^c)). This de- 
fines G?, G^, G^ and G^, their associated photoelastic 
signal [2g. We then estimate the normal force of each 
link, Fjv, by Fn = {Gl + Gl + Gl+ G^)/2. In the same 
vein, we estimate the tangential force of each link. Ft, 
by Ft - {Gl ~ Gl + Gl - Gl)/2. 



C. Calibration and units 

We compute G^, the average of Gf over space, and 
compare it with the force F, measured by the force sen- 
sor (d), normalized by Mg, the total weight of the grain 
assembly (figuredja)). One observes a linear relationship 
between the force measured at the piston and the sum of 
Gj over the entire picture, G^. In the following, we use 
this same linear relation to calibrate the local Gf . Below, 
all pressures and forces computed using the photo-clastic 
images, arc expressed in units of Mg. Lengths are ex- 
pressed in units of the small grain diameter and time is 
expressed in units of the microscopic time determined by 
the stiffness of two compressed discs: to = (fc/m)~^/^, 
where m is the mass of a grain (~ 3.75 x 10~^ kg) and k 
is the stiffness of two compressed disks ('^ 1.5 x 10"^ N/m). 



D. Protocol: Obtaining a granular glass 

As already emphasized in the introduction, the Jam- 
ming transition is intrinsically a T = 0, out-of- 
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equilibrium transition, and therefore depends on the pro- 
tocol followed to prepare the packing of interest. The 
situation need not be made simpler by the introduc- 
tion of thermal or mechanical vibration. Indeed, for 
the packing fractions of interest, most systems become 
naturally dynamically arrested in non-equilibrium glassy 
states. The steep increase of the relaxation times associ- 
ated with glassy behavior seriously hampers experimen- 
tal work [29l - l32| : samples brought to the high packing 
fractions of Jamming are deep into the glass phase and 
are difficult to manipulate on reasonable timescales. For 
athermal granular media, the situation is similar: they 
need some mechanical energy to be maintained in a non- 
equilibrium steady-state (NESS). As for thermal systems, 
this requires extremely slow compaction of the sample 
in order to avoid aging dynamics on the experimental 
timescales [s^ H^l- For that reason, most granular ex- 
periments actually probe the glass transition and not the 
Jamming transition [ssHstI . 

Here, we perform an annealed compaction (figure |4}d) , 
i.e., we increase packing fraction by constant amounts of 
(5(^ = 3 X 10"'*, with exponentially increasing time steps. 
Then, the packing fraction is stepwise decreased, and 
measurements are performed between the decompaction 
steps (figure dfa-b)). Lechenault et al. Q checked that 
the dynamics is reversible and stationary on experimen- 
tal time scales during these decompression steps. As sug- 
gested by figure [3 the structure of the packing we obtain 
following the above protocol is frozen: the superimpo- 
sition of two Voronoi tessellations, separated by a time 
lag of 5000 vibration cycles, display very few rearrange- 
ments, even for the lowest packing fraction. Such rear- 
rangements are further quantified by Q™{t), the aver- 
age fraction of neighbor relationships surviving in a time 
interval t. Plotted with respect to the lag time, r, Q™ 




FIG. 5: Obtaining a granular glass, (a): Superposition 
of the Voronoi' cells computed at times t = 1 and t = 5000 
for the loosest packing (0 = 0.8031). (b): Average fraction 
of neighbors Q"'"'{t) which have not changed between two 
images separated by a time interval r, for different packing 
fractions. The vibration frequency is 10 Hz and the packing 
fraction takes 13 values in the range [0.80 — 0.82]. The color 
code spans from blue (low packing fractions) to red (high 
packing fractions). 



remains larger than 95% even for the loosest packing frac- 
tion, and barely departs from 1 for the densest ones {Hg- 
ure[5l^b)). In the language of the glass community, "there 
is no a relaxation" , meaning that the density profile sur- 
vives on the experimental time-scale and the system can 
safely be considered as a glass, the structure of which is 
essentially frozen. 

Finally, note that despite the fact that we perform the 
same protocol for each experiment, the initial conditions 
are still different for each run. Also, the system size is fi- 
nite, and therefore, the Jamming transition of each pack- 
ing will fluctuate from one realization to another. It is 
important to keep this in mind when comparing indepen- 
dent experimental runs. 



III. PRESSURE AND CONTACT FORCES 

For the ideal case of soft spheres at zero temperature, 
the pressure inside the packing exhibits the same basic 
features as the energy: below Jamming, it is strictly zero 
and above Jamming it grows with the packing fraction, 
according to the interaction force between particles. It 
is thus of interest, as a first sight at the transition in a 
system with dynamics, to look at the dependence of the 
pressure with the packing fraction. 

Figure [6] displays the pressure measured at the wall as 
a function of the packing fraction. Ptot (respectively 
Pstat) is the pressure measured when the vibration is 
applied (vibration on) or not (vibration off). Pstat cor- 
responds to the the static pressure sustained by the pack- 
ing whereas Pdy n = Ptot — Pstat is the dynamic part 
of the pressure that comes from the vibration. One ob- 
serves a smooth crossover from a constant, but nonzero 
pressure, to an pressure that increases with the pack- 
ing fraction. On the large packing fraction side of the 
crossover, Ptot — Pstat and the pressure, which is 
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FIG. 6: Wall pressure vs. packing fraction: (Q) : Ptot, 
(□): Pstat, (A): Pdyn , as defined in the text for the present 
PSM-4 disks experiment. Note: the finite stiffness of the 
piston has been calibrated and removed from the data. The 
vibration frequency is / = 10 Hz. 
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mostly static, follows what is expected from the zero tem- 
perature prediction: it increases with packing fraction, 
according to the particle stiffness. On the low pack- 
ing fraction side, there is an irreducible kinetic part of 
the pressure, induced by the vibration. To zeroth order, 
the crossover corresponding to Jamming can be identi- 
fied with the packing fraction where the static pressure 
becomes larger than the kinetic one. Note that the static 
part of the pressure is not strictly zero below the cross- 
over. We attribute this to the mobilization of the friction 
at the contacts, when the vibrating board is stopped. We 
return to the possible roles of friction in the discussion 
section. 

One must realize that the kinetic part of the pres- 
sure, which is observed on the loose side of the Jamming 
crossover, does not strictly speaking come from collisions 
of the grains with the wall. Indeed, the instantaneous 
dynamics is very different from that of a thermal liquid, 
where the pressure has a coUisional origin. Here, the forc- 
ing is periodic and a priori strongly anisotropic. The par- 
ticles are accelerated along the vibration axis, then com- 
pressed along one wall, before being accelerated back in 
the reverse direction. A clearer idea of this process comes 
from the dynamics of the average inter-particles forces 
during a few vibration cycles, and by a decomposition 
of these forces in to the vibrat ion and tran vserverse di- 
rections: Fx = \/ {{fij ■ exf) and Fy = \J {{fij ■ eV)^), 
where (•) is the average over space, and where ex and eV 
are unit vectors along the vibration and the transverse 
directions. 

One observes in figure Elja) that for the low packing 
fractions, there are strong oscillations at the vibration 
frequency. These oscillations correspond to the compres- 
sion of the grains on the side walls. Interestingly, these 
oscillations are in phase, within the temporal resolution 
of the acquisition: the transfer of momentum, from the 
direction of vibration to the transverse direction, is in- 
stantaneous, as compared to the time scales considered 
here. 




(a) (b) 

FIG. 7: Short time photo-elastic response, (a): Av- 
erage interparticle force decomposed into the vibration Fx 
and the transverse Fy directions vs. time t for packing frac- 
tions (j) = 0.8079 (blue), 0.8123 (green) and_a_8196 (red), 
(b): Time averaged quantities: Fx (blue), Fy (red) and 
F — (Fx + -FV)/\/2 (green) vs. packing fraction 4>. The 
vibration frequency is / = 10 Hz. 



For the packing fraction above the kinetic to static 
crossover, the oscillations are much less pronounced, and 
the periodicity not so clearly defined: in that regime, 
the global motion of the grains with respect to the os- 
cillating plate is reduced. Finally, averaging temporally 
those signals, and plotting the averages as functions of 
the packing fraction (see figure [Tl^b)), the same trends 
occur as those observed for the pressure measured at the 
wall, albeit with larger fluctuations, since the temporal 
sampling is much smaller. Note that the average force in 
the direction of vibration is only slightly larger than the 
average force in the transverse direction, indicating that 
the redistribution of momentum ensures the formation 
of a rather isotropic force network. We confirm this by 
a direct inspection of the pressure field inside each grain 
Gi, interpolated on a cartesian grid, and the computation 
of the spatial auto-correlation function (figure [U . They 
both confirm a good level of isotropy of the pressure dis- 
tribution in the packing. 

We conclude this section by inspecting the proba- 
bility distribution of the normal and tangential inter- 
particle forces Fpf and Ft defined by fij = F^fij + 
Fxtij, where fy (respectively tij) is the normal (re- 
spectively tangential) vector between the two grains 
i and j. To be precise, we consider the distribu- 
tion of = [Fjv(i,t)/(G(i,t)),](G(i,t)),,t and Ft = 
[FT{i,t)/{G{i,t))i]{G{i,t))i^t, where < . >i is the instan- 
taneous average over the particles and < . >i^t is the 
avergae over time and space. Such a normalization has 
the advantage of capturing the shape and width of the 
distributions, without including the temporal variability 
of the packings Q, hence avoiding spuriously large tail 
distributions. Also, since the normalization is the same 
for Fm and Ft, the ratio of Ft/Fn = Ft/F/N, which 
ensures a correct interpretation in terms of friction co- 
efficient. The distributions (figure ^a) and (b)) have 
exponential tails at all packing fractions, and widen as 
the packing fraction is increased. This is consistent with 
existing works on granular packings [H, HI]. Note, how- 
ever, that the existing consensus on the exponential tails 
of force distributions is not founded on any unambiguous 




(a) (b) 

FIG. 8: Isotropy of the force network, (a): Interpolated 
instantaneous Gi's on a cartesian grid and (b): its associated 
2d-autocorrelation. The packing fraction is (j> ~ 0.82. The 
vibration frequency is / = 10 Hz. 
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FIG. 9: Distributions of forces, (color online) Distribution 
of normal (a) and tangential (b) forces. The vibration fre- 
quency is / = 10 Hz and the packing fractions are the same as 
in figure [5] Joined distributions of tangential Ft vs. normal 
Fn forces for (c): (j) = 0.8178, and (d): (f) = 0.8125; same 
vibration frequency. 

arRuments, and that some studies report non-exponential 
tails [23, Sll . Leaving aside this debate, we choose 
to focus on the joint distributions of Fm and Ft {^^g- 
urc[5]Jc) and (d)). The ratio Ft/Fn is close to 0.2. on 
average, and always smaller than 0.7, which provides a 
good estimate for the static friction between the PSM- 
4 disks. One also notices an accumulation of contacts 
close to the threshold value fig, especially at low forces, 
where a gap in the distribution clearly separates a major- 
ity of contacts with Ft/Fn — 0.2 from a secondary peak 
of contacts with Ft/Fm ^ Ms — O-'''- These so-called 
"critical contacts" arc on the verge of slipping. Whether 
these slipping events are trivial fluctuations, or contain 
some interesting correlations in the vicinity of the Jam- 
ming crossover, was the central issue discussed recently 
by the present authors [2^. This issue will be recast in 
section HVBl 



IV. DYNAMICS OF THE CONTACT 
NETWORK 

In order to measure 2i(t), the number of contacts of 
particle z at a given time i, one must identify the potential 
contacts of particle i with its neighbors by thresholding 
the normal force, fjy and the inter-particle distance, s. 
We have shown [26| that the overall behavior, including 
the statics and the dynamics, remains unchanged when 



varying the thresholds within a reasonable range. Here, 
we shall keep the threshold fixed, and focus on the dy- 
namics of the contact network both at short and long 
times. 



A. Statics and short time dynamics 

Figure fTUT a) shows the average number of contacts z, 
computed over the images acquired in a stroboscopic way 
at a vibration frequency of 10 Hz, vs. the packing frac- 
tion. These data exhibit a clear cusp at any given packing 
fraction, (f)^ . For packing fractions larger than z in- 
creases with the packing fraction, in a way that is similar 
to what is reported for zero-temperature soft spheres. By 
contrast, for packing fractions smaller than z \s non- 
zero constant, and there is no discontinuity across (f)\ 
in contrast with the zero temperature behavior of soft 
spheres. 

In terms of the pressure signal, a more precise pic- 
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FIG. 10: Static and short time dynamics of the con- 
tact network (a): Average contact number obtained from 
the stroboscopic data, vs. packing fraction 0. The dashed 
line indicates (j>^ = 0.8151. (b): Instantaneous average con- 
tact number z vs. time t, at (p = 0.8079 (blue), 0.8123 (green) 
and 0.8196 (red). The green dashed lines indicate the times, 
tl;, where the contact number is minimal, namely when the 
grains experience the smallest acceleration. The red dotted 
lines indicate time frames t| where, by contrast, the grains are 
compressed against a wall, (b): Different temporal averages 
of the contact number as a function of the packing fraction: 
in blue is the contact number averaged over all acquisition 
frames; in green, respectively in red, is the contact number 
averaged over the time frames tl, respectively t|. The vibra- 
tion frequency is / = 10 Hz 
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ture of the mechanisms at play behind the shape of the 
z{4>) dependence can be obtained by examining the dy- 
namics during a vibration cycle. Figure fTOTb) displays 
the instantaneous contact number, z{t) = 
acquired with the fast camera for three different pack- 
ing fractions. Here, N is the number of particles. For 
lower packing fractions, strong oscillations at the vibra- 
tion frequency are clearly visible, while they are reduced 
and not so well defined at larger packing fractions. The 
similarity with the force signals reported in figure is 
striking, and one easily understands that the number of 
contacts is temporarily larger when the grains are com- 
pressed against the wall. As a result, the average number 
of contacts computed from the stroboscopic data depends 
on the precise phase at which the acquisition is performed 
and this dependence is most significant when the packing 
fraction is low. This is illustrated in figure rTOf c). where 
temporal averages of the contacts number, acquired at 
different phases, are plotted as a function of the pack- 
ing fraction. In green, is the contacts number averaged 
over time frames which arc in phase with the minimal 
acceleration: the grains are "away" from the walls, and 
the number of contact is minimal too. This situtation 
corresponds to the vertical green dotted lines in figure [TU] 
(b). There is no longer any evidence of a cusp, and the 
crossover is only indicated by an inflexion point that is 
barely discernable. In red, is the contact number aver- 
aged over time frames where the grains are compressed 
on one of the walls and the number of contacts is maxi- 
mal; this situation corresponds to the vertical red dotted 
lines in figure [TUf b). One recovers the cusp observed in 
figure llOf a) , for which the stroboscopic acquisition was 
indeed performed in phase with the maximal accelera- 
tion and the minimal velocity of the plate. As explained 
in section |lll this choice of phase minimizes blur in the 
images 

Altogether, the presence of the cusp is related to the 
vibrational forcing and to the specific phase at which the 
stroboscopic acquisition is performed. There is no singu- 
larity in the dependence of the average contact number 
as a function of the packing fraction. As for the pres- 
sure, the signature of the Jamming transition is replaced 
by a crossover, the precise location of which depends on 
the details of the measure. Here, we were lucky enough 
to capture the images at the phase of the vibration, for 
which the crossover t/)^ is most easily identified (red curve 
of figure [TUb). One should however keep in mind that the 
green curve in the same figure is actually a more realistic 
dependence of the average number of contacts with the 
packing fraction [5|. 



B. Long times dynamics 

The results of this section, which concern the long 
time dynamics of the contact network and the nontrivial 
correlations that it contains, have been reported previ- 
ously (26j . Here, we recast the important message they 
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FIG. 11: relaxation dynamics of the contact network. 

(color online), (a): Temporal average of the overlap function 
Q^{t). Packing fractions as in figure O (b): Plateau value 
of Qs, defined by = Q^{t = 4) (a), and fraction of non- 
rattling particles ( X ), vs. packing fraction (/). (c): Relaxation 
time of the contact network, (see text for definition) (■, 
left axis), and average contact number, z {x, right axis), vs. 
packing fraction, (j>- The plain red line is a fit of the form 



At 



The dashed line indicates 



At 



0.8151. 



The vibration frequency / = 10 Hz. 



convey, together with new results. 

To characterize the dynamics of the contact network, 
we introduce the contact overlap function, which evalu- 
ates how much the contacts have fluctuated between t 
and t + t: 



Q'{t,r) = -J2Qt{t,r), 



i=l 



where 



1 if |z,(t + r) -z,(t)| < 1 
Oif |z,(t + r)-z,(t)| >1 



(1) 



(2) 



Other choices of overlap functions are possible, and have 
been tried: the present results do not crucially depend on 
the particular choice. Figure ITlT a) displays the temporal 
average of Q^{t, r) for a vibration frequency of 10 Hz and 
the same set of packing fraction (f> G [0.80 — 0.82] as in the 
above sections. At rather large packing fractions, (5^(t) 
is constant, with a plateau value, which depends weakly 
on the packing fraction. Hence, there is no relaxation 
on long time scales of the contact network. The relax- 
ation, which occurs at short times, and is responsible 
for the plateau value, cannot be observed in the present 
stroboscopic data. However, it is apparently related to 
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the motion of the ratthng particles, i.e., particles hav- 
ing less than 2 contacts, as suggested by the very strong 
correlation observed between the fraction of non-rattling 
particles and the value of the plateau at short-times (see 
figure [TlT b)). At lower packing packing fractions, a long 
time decorrelation sets in. We define the relaxation time 
of the contact network, r^, such that Q^(t^) = 0.9. Note 
that this value of 0.9 is rather large as compared to the 
most commonly used value of 0.5. However, it is the 
smallest one which allows the measure of for a broad 
range of packing fractions. We note that relaxation times 
measured in a standard way would be orders of magni- 
tude larger. As shown in figure [TTJc), left axis, in- 
creases sharply with the packing fraction, and possibly 
diverges at the packing fraction (ji^ = 0.8151, where the 
average number of contact starts to increase with the 
packing fraction. 

Interestingly, the dynamics of the contact network be- 
low 0^ exhibits strong fiuctuations and dynamical hetero- 
geneities, albeit of a different kind from those reported 
in the literature, when studying the dynamics of super- 
cooled liquids close to their glass transition (see [lal)- 
Here, the heterogeneities are relative to the degrees of 
freedom describing the contacts, not the position of the 
particles. To quantify such heterogeneities, one can com- 
pute the dynamical susceptibility which estimates the 
range of spatial correlations in the dynamics of the con- 
tact network: 



Var{Q^{t,T)) 
{Var{Ql{t,T))),' 



(3) 



where Var[.) denotes the variances sampled over time 
and denotes the average over the grains. Y|(r) has a 
maximum for r = r* (not shown here, see [26[), and we 
have studied how the maximum x|* of Xii'^) depends on 
both the packing fraction and the vibration frequency. 
To do so, it was necessary to prepare different packings, 
and run independent experiments at three different vi- 
bration frequencies: / = 6.25, 7.5 and 10 Hz. As em- 
phasized in section ini the precise value of the Jamming 
transition, and certainly those of the crossovers reported 
here, depend on the specific packing. Hence, following 
the methodology of section IIV Al for each frequency: we 
identified the structural crossover 0^ , from which we de- 
fine a reduced packing fraction e = (0 — </>''')/(/)''', in order 
to compare the different experimental runs. Note that 
for frequencies smaller than /g = 4.17 Hz, the grains 
do not slip on the driving plate, and the mechanical 
excitation is effectively null. Accordingly, we introduce 
7 = (/ - /o)//o, 7 = 0.5, 0.8, and 1.4, to quantify the 
level of mechanical excitation. 

The results are summarized in figure [T^] and can be 
found in more detail in [2^. x\* is non-monotonic with 
respect to the reduced packing fraction, and has a max- 
imum value at a negative reduced packing fraction e*. 
This indicates the existence of a dynamical crossover cor- 
responding to a maximally collective relaxation of the 
contact network at a packing fraction lower than the 




FIG. 12: Towards zero vibrations, (color online). Maxi- 
mal dynamic susceptibility of the contacts, X4*i vs. reduced 
packing fraction, e and reduced vibration magnitude, 7. 



structural crossover. Also, when 7 is decreased one ob- 
serves that («) e* vanishes, i.e., the location of the dy- 
namical crossover moves towards and {ii) the mag- 
nitude of the maximum Xa* significantly increases as 
1/7. Hence, we can safely conjecture that in the limit 
of no effective mechanical excitation the structural and 
dynamical crossovers merge, while the length scale as- 
sociated with the dynamical crossover diverges. This 
strongly suggests, that we have probed the vicinity of 
a critical point, which in the present case ought to be 
the Jamming transition in the absence of dynamics. As 
a matter of fact, a similar phenomenology occurs for 
equilibrium systems close to a thermodynamic critical 
point: at the critical point, thermodynamic susceptibil- 
ities diverge and away from it, in the supercritical re- 
gion, they exhibit finite maxima. These are the so-called 
Widom lines [i^, . Recent experiments [13, prob- 
ing the dispersion of nano-metric acoustic waves report 
a crossover of the acoustic dispersion along one of such 
Widom lines and demonstrate the existence of a dynami- 
cal crossover involving subtle mechanisms at the particle 
scale in the supercritical region a thermodynamical crit- 
ical point [46j . 

The above results clearly indicate that the mechani- 
cal agitation blurs the singular nature of the Jamming 
transition. This is a similar effect to one reported in the 
presence of thermal agitation for the Jamming transition 
of soft spheres [l| . One of the remarkable results of that 
work is that the authors demonstrate in a convincing 
manner that all the physics of the soft sphere systems 
close to Jamming can be captured by a careful examina- 
tion of the mean square displacement of the particles as 
a function of time. They can then use this measure as 
an effective thermometer to locate existing colloidal ex- 
periments in the Temperature-Packing fraction parame- 
ter space. Is it possible to extend the approach to the 
present case of vibrated granular media? If yes, where 
does the present experiment sits in an equivalent param- 
eter space? 
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(a) (b) 

FIG. 13: Motion of the center of mass, (color online) 
(a) Center of mass position, in the vibration direction Xi, 
(top), and in the transverse direction Yt (bottom), vs. time 
t,at(j} = 0.8089 (blue), 0.8161 (green) and 0.8196 (red), (b): 
Amplitudes (+) and (O) vs. packing fraction, ij). The 
vibration frequency / = 10 Hz. 



V. DISPLACEMENTS FIELDS 

In order to answer the above questions, one needs to 
extract the mean square displacement of the particles 
on the largest possible range of timescales. While this 
is a straightforward but CPU costly task in numerical 
simulations, we shall see that it requires rather intricate 
data analysis in the present experiment. The reasons are 
twofold. First, the short time and long dynamics arc ac- 
quired independently and in very different ways. While 
the long time dynamics is acquired in phase with the vi- 
bration, the short time dynamics are acquired within a 
vibration cycle. The long time acquisition naturally fil- 
ters the "trivial" motion of the plate, but the short time 
does not, and we will have to filter it out. Second, we 
shall see that on long time scales, low-amplitude convec- 
tion occurs. Although the resulting flow is never large-it 
mostly consists of a non-monotonic solid body rotation- 
we shall remove it before computing the mean square 
displacement. 



A. Short time oscillations 

The motion of the center of mass {Xb{t),Yi,{t)) ^ 
{{Xi{t))i, {Yi{t))i) provides a good indication of the way 
the energy is injected in the system at large scale. Fig- 
uresfTSTa) shows that the center of mass oscillates period- 
ically, with a period equal to the forcing frequency. The 
amplitude of the motion is much larger in the direction 
of vibration, but part of the forcing is transfcred to the 
transverse direction too. As shown in figure [THTb). the 
amplitudes of the oscillations, and A^ depend on 
the packing fraction: they are constant at low packing 
fractions, typically when < 0*, and they decrease for 
larger packing fractions, suggesting that energy injection 
is less efficient at large packing fractions. 

In order to investigate the way the energy is trans- 
fcred to smaller scales, we compute the averaged spectral 
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(a) (b) 
FIG. 14: Spectral properties, (color online) (a): Average 
Fourier Energy Spectral Density Xj (blue) and Y^ (red) of 
the grain position fluctuations for (f> = 0.8196. The dotted- 
dashed lines indicate the first ten harmonics of the excitation 
frequency, (b): Average Fourier Energy Spectrum Density 
= Xf + Yf of grain positions fluctuations at cj> = 0.8089 
(blue), 0.8161 (green) and 0.8196 (maroon), after filtering the 
trajectories as described in the text. Inset : Low frequency 
limit (/() = 0.3), of the Average Fourier Energy Spectrum 
Density, A^j, vs. reduced packing fraction, e = {(fi — (j>^)/cj>^ . 
The vibration frequency / = 10 Hz, i.e 7 — 1.4. 



density of the positions fluctuations. Specifically, we de- 
fine (X.it) = X,{t) - Xb{t), Y,{t) = Y,it) - Ybit)), corre- 
sponding to the grain trajectories in the frame of refer- 
ence of the oscillating center of mass. We next compute 

Xi'if) = ESD{X,{t) - {X,it))t), and similarly f,'(/), 
where ESD denotes the Fourier energy spectral density 

(ESD). We then average over the grains to obtain the 

^ 2 ^ ^2 

spectra Xj = {Xi {f))i and Y^ = {Y^ (/)),;, displayed in 

figure [Ula). 

The energy cascades down to high frequencies, which 
unfortunately preserves the signature of the periodic forc- 
ing in the form of strong harmonics. This indicates that 
considering the motion of the grains in the frame of the 
center of mass is not sufficient to completely filter out 
the periodic motion induced by the moving plate. We 
thus further filter the grain trajectories by applying a 
band-cut Butterworth filter centered on each harmonic 
(up to the fifth) and a low-pass Butterworth filter with a 
cut-off frequency of 5 times the vibration frequency, on 
Xi{t) and Yi{t). The spectra of the filtered trajectories. 
Ay- = Xj + Xj (figure lMT b)) confirm that the harmonics 
have been successfully filtered out. The resulting motion 
A^° at the lowest frequency (/o = 0.3), corresponding to 
a timescale of a few vibration cycles, is a good indica- 
tor of the typical cage size in which the particle vibrates. 
It strongly depends on the packing fraction and sharply 
decreases as the Jamming crossovers are crossed. The 
absolute magnitude of (10^^ to 10^"*) corresponds 
to a typical cage size of ^ 10~^ to 10~^ grain diameter. 

In the remainder of the paper, we will apply the fil- 
tering procedure described here on each grain trajectory 
prior to computing any statistical property for the fast 
camera data. 
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FIG. 15: Solid body rotation, (color online) (a): Grain 
displacements over a lag time r = 6 x 10*^ at e = -0.0948 
(inset: zoom on a few grains at the edge of the region of 
interest), (b): Orthoradial displacements vs. distance to 
center for a lag time r = 6 x 10® at e = —0.0948, and 7 = 1.4. 
(c): 51/r vs. time, t, for different lag times r = lO'' (blue), 
T — 10" (green) and r = 10^ (red) at a reduced packing 
fraction, e = —0.0948. (c): Rotational drift coefHcient SI vs. 
reduced packing fraction, e. The vibration frequency / — 
10 Hz, i.e 7 = 1.4. 



B. Long time rotation 

We now turn to the stroboscopic trajectories of the 
grains. Figure IWa) displays the displacement of all 
grains in the region of interest (ROI), integrated over 
a lag time t = 6 x 10®. The inset provides a zoom 
on the trajectories of a few grains at the edge of the 
ROI. One observes a clear global rotation, which, curi- 
ously and fortunately, is essentially solid body motion, 
as demonstrated by the linear dependence of the ortho- 
radial displacement Ri{t,t + T){9i[t + t) — 9i{t)) with 
the distance Ri{t,t + r) = |[fi(£±2±tnl£)|| to the center 
of the cell (figure [TSlb)). It is fairly easy to remove 
this solid body rotation from the grain displacements 

Kn{t) = n{t + T)-nit). 

One defines : 



where 
A?rl(0 



n{t)T\ / r-.(t)-Hn(t+r) ^ 

-n{t)T J \ 2 'r^tj 



(4) 



(5) 



angular velocity J7(i) and the center of rotation r^{t) are 
explicitly computed from the displacements AT-r,;(t), by 
minimizing {\\fi{t + t) - ri{t) - Ar, iTt{t)\\'^)i, with re- 
spect to il{t) and r^{t). One finds: 



r,{t)+ri{t+T) 
2 



and 



(6) 



N / \ 

r?(<) = E ( -1 0) (^^'(^ + ^) - W) • 



Figure fTSl c) reveals that x T{t), the angular rotation 
between times t and t + r, fluctuates around zero, mean- 
ing that the solid body rotation has no prefered direction. 
As a result, there is no statistically systematic drift in any 
direction. However, for any finite time interval, [<,t-|-T], 
there is a finite angular displacement, the magnitude of 
which is controlled by As shown in Figure [TSl d). it 

sharply decreases as the packing fraction increases across 
the Jamming crossovers. 



C. Resulting vibrating dynamics 

Now that both the short time "trivial" oscillating mo- 
tion and the long time convection have been filtered out, 
we are in a position to characterize the vibrating dynam- 
ics of the grains in the frozen structure of the packing on 
time scales ranging from a hundredth of a cycle to several 
thousands cycles. We compute the following estimator of 
the mean square displacement: 



MSD = - {{\Arr 



(8) 



is the solid rotation deformation field. The values of the 



where (...) denotes the average over times and parti- 
cles, and Arr is the particle displacement obtained from 
the filtering procedures described in the previous section. 
The choice of this estimator is motivated by the fact that 
it ensures a lower statistical weight to very large moves, 
such as those of the rattling particles. The factor ^ en- 
sures quantitative matching with the proper mean square 
displacement in the case of gaussian statistics. Alterna- 
tively, one could remove the rattling particles, but that 
strategy requires additional filtering and or thresholding. 
Figure [in] displays the mean square displacement over the 
full timescale interval probed in this experimental study. 
We again emphasize that the data at short times, shorter 
than 10'^, were obtained from the fast recording of the 
grain motion within the vibrating cycles, while those at 
long times were obtained performing stroboscopic acqui- 
sition in phase with the oscillating driving plate. Each 
type of acquisition were performed during independent 
experimental runs. The color codes the packing fraction. 
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(b) (c) 
FIG. 16: Mean square displacements, (color online), (a): 
Mean square displacements MSD for filtered trajectories (see 
text) vs. lag time r for both the short time (fast camera) 
and long time (stroboscopic acquisition) experiments. The 
packing fractions explore the same range and are color coded 
as in figure [5] the binning being finer for the fast camera 
acquisition, (b): Plateau value obtained from the short 
time data MSD (+), from the long time stroboscopic data 
MSD ( X ) , and from the low frequency limit of Energy Spectral 
Density, A^^ (□) vs. reduced packing fraction, e — {(j> — 
cj>^)/4>^. (c): Plateau entrance time rsn vs. reduced packing 
fraction e = {(fi — (f)'' ) / (f>^ . The vibration frequency / = 10 Hz, 
i.e 7 — 1.4. 



The good overlap of the mean square displacement at in- 
termediate time scales is not enforced and is remarkably 
good. 

Altogether, one observes three regimes: a ballistic 
regime at short time r < teu, a plateau at intermediate 
time scales, teu < t, and for low enough packing fraction 
a crossover towards a diffusive regime at long time scales. 
The plateau regime characterizes the vibrational dynam- 
ics we are interested in. The height of the plateau, A^, 
measures the square of the average vibration amplitude 
of the grains within their cage. It decreases from 10""' to 
10^^ for increasing packing fractions (figure HBT b)). and 
it is consistent with the first estimate of the cage size, 
we had obtained in section FV Al from the low frequency 
limit of the Fourier spectral density of the position fluc- 



tuations, A^y. The short time entrance to the plateau, 
estimated by teu = (A^/ivT)!/^, where K ~ 10"^ is ob- 
tained from the analysis of the ballistic regime, typically 



occurs at teu 100 and slightly decreases as the packing 
fraction is increased (figure [W c)-left axis); the larger the 
packing fraction, the sooner the grains feel their neigh- 
bors and enter the vibrational regime. 

The above vibrational dynamics is very similar to the 
one reported for thermal harmonic sphere systems close 
to Jamming [l|, \^ . In this later study, a ballistic regime 
occures at short time, followed by a plateau regime, the 
height of which decreases strongly with the packing frac- 
tion when crossing the Jamming point. A plateau exit 
is also reported in , where the authors show that the 
plateau exit time increases when the quench rate used to 
prepare the packing is decreased. This plateau exit is not 
reported in [l|. However, the maximum lag time there 
was 10^, and the systems were carefully equilibrated, so 
that the plateau exit, if it existed, was probably much 
larger than the simulated timescales. Before addressing 
the more quantitative comparison, which will allow us to 
discuss whether thermal soft spheres are a good model 
for mechanically excited grains, we will finish the descrip- 
tion of the dynamics by characterizing its heterogeneities. 
Note that these heterogeneities, first reported in the brass 
grains experiment [23| and more recently in the harmonic 
spheres simulation [l| are distinct from those encountered 
in super-cooled liquids when approaching the glass tran- 
sition [i^. Here, the structure is frozen, hence, the het- 
erogeneities are not related to the relaxation of the struc- 
ture. The next section will show how they are related to 
the heterogeneities of the contact dynamics described in 
section HFEl 



VI. DYNAMICAL HETEROGENEITIES 

In this section, we investigate the heterogeneities of the 
particle displacements. To do so, we focus on the long 
time stroboscopic data, once the convective motion has 
been subtracted. We will show that these heterogeneities 
take place at very small scales and are temporally cor- 
related to the heterogeneities of the contact dynamics. 
Finally, a closer look at the organization of the contacts 
at short time will demonstrate that these heterogeneities 
take their root in the short time organization of the con- 
tact network, namely in the vibrational dynamics of the 
structure. 



A. Heterogeneous non-afflne dynamics 

The characterization of dynamical heterogeneities has 
now become a standard tool in the study of the dynam- 
ical slowing down of super-cooled liquids and/or colloids 
approaching their glass transition [48[ . It is much less fre- 
quently used when probing the Jamming transition, but 
relics on the same procedure fl6j . In order to character- 
ize the dynamics, and in particular to probe collective 
effects, one defines a dynamical structure factor for the 
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(c) (d) 
FIG. 17: Dynamical heterogeneities, (color online) (a): 
Dynamic susceptibility of the displacements xK^**! ''') "^s. lag 
time r. Same packing fractions as in figure (5] (b): Maximal 
dynamical susceptibility of the displacements xl vs. reduced 
packing fraction e. (c): vs. reduced packing fraction e. 
(d): a* vs. reduced packing fraction e. Dashed lines in frames 
(b,c,d) indicate e*. The vibration frequency / = 10 Hz, i.e 
7 = 1.4. 



displacements, AT-ri(t) (defined in eq. 2]): 
Q^lt,r,a) = l^Qf(t,T,a), 



(9) 



where 



(10) 



Qf(t,T,a) = exp(-||A,rl(t)||V2a'). 
This dynamical structure factor probes the dynamics 
at scale, a, and time t: Ql{t,T,a) ~ 0(1), when the 
particle i has moved more (less) than a, during r. One 
then defines the dynamic susceptibility: 



N 



Var{Q^{t,T,a)), 
(11) 

where Var denotes the temporal variance. It provides an 
estimate of the average number of particles which move 
up to the distance a during a time r in a correlated man- 
ner. In general Xii'^j^) absolute maximum 
for a = a* and t = (see for instance j23j). 

Figure [T7] illustrates the way the dynamical hetero- 
geneities depend on the packing fraction. The most im- 
portant effect is that is nonmonotonous and exhibits 
a clear maximum at precisely the reduced packing frac- 
tion e* (figure [TTT b)). The magnitude of X4 close to 
e* is close to 100, roughly a tenth of the total number 



of particles. (Even closer to e*, one may note the data 
point indicated by □, which is anomalously low com- 
pared to the trend given by the other data points. We 
believe that this is a signature of the lack of statistics 
necessary to resolve much larger heterogeneities.) The 
timescale, , where this maximum occurs, is not very 
sharply defined (note the logarithmic scale for r), as 
can be seen from the depcndancc of Xiio,*^''') on r (fig- 
ure [TWa)). But it clearly increases significantly when 
the packing fraction increases and certainly is larger than 
the timescalcs for packing fractions larger than 0^ (fig- 
ure [T71[c)). The length scale, a*, over which the parti- 
cles move while building up the heterogeneities, decreases 
with the packing fraction and is of the order of 10~^(i (fig- 
ure [TTJd)). The same observation made in the case of the 
brass disks [1^ lead the authors to conclude that the dy- 
namical heterogeneities observed close to Jamming have 
their origin in the dynamics of the contacts. We are now 
in position to confirm this intuition. 



B. Relation to contact dynamics 

The fact that the heterogeneities observed in the dy- 
namics of the contact and in the displacement field are 
both maximal at the same value of the reduced pack- 
ing fraction e* is already a strong indication that they 
have a common origin. This is further confirmed by the 
quantitative comparisons of Xi lo xl* ^-nd of to 
T^* provided in figure [TST a) and (b). Xi ^-nd are 
respectively proportional to xl* ^-iid t^* confirming a 
strong correlation between the two aspects of the dy- 
namics. As in figure [TW b). the □ data point is way off 
the trend given by the other data points, because of a 
lack of statistics at the Jamming crossover (see previous 
section). Whereas the timescales are essentially identi- 
cal, the dynamical susceptibility associated with the dis- 
placements is 20 times larger than that associated with 
the contacts. One must remain cautious in the inter- 
pretation of such a factor, since the dynamical suscepti- 
bilities are only an indicator of the number of elements 
correlated, even when they are properly normalized by 
the intrinsic fluctuations: the shape of the spatial cor- 
relator also enters into play. With that caveat, such a 
large factor suggests that the spatial organization of the 
dynamics is different in the two cases. This is indeed con- 
firmed in figureflSTc) and (d), which shows the snapshots 
of respectively Q^{t,T*) and Q\{t,T* ,a*), taken at the 
same time. Whereas the dynamical heterogeneities of the 
displacements are organized in well identified large clus- 
ters, those of the contacts seem more scattered in smaller 
chain-like clusters. The dynamical correspondance is not 
simply that the particles moving more than a* lose or 
gain contacts. On the contrary, it suggests that the loss 
of contact at some place induces motions on the scale of 
a* further away, and in turn, the loss of other contacts. 

In section HVBl we have seen that dynamical hetero- 
geneities of the contacts enlarge when the vibration frc- 
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(a) (b) 




(e) (f) 

FIG. 18: Dynamical heterogeneities, (a): Maximal dy- 
namical susceptibility of the displacements xJ vs. maximal 
dynamical susceptibility of the contacts xl* a^nd (b): vs. 
T^* in parametric plots, where each point correspond to a dif- 
ferent packing fraction (same color code as in figure [5]). The 
vibration frequency / = 10 Hz, i.e 7 = 1.4. (c-d): Maps 
of gf(t,r) (c) and Qf(t,r,ai/2) (d), for e = -0.0013. (d): 
Color code spans from yellow (Qf {t,T) = 0) to red (Qf (t,T) = 
1). (d): Color code spans from blue (QI (t, t, ai/2) = 0) to 
red (Qi (t, T, ai/2) = 1). The vibration frequency / = 10 Hz, 
i.e 7 = 1.4. (e): Same plot as in (a) but for different values 
of 7 = 0.5 (blue), 0.8 (green), and 1.4 (red), (f): Peak of the 
dynamical susceptibility of the displacements max(x4 ) and 
peak of dynamical susceptibility of the contacts max(x4*) vs. 
7- 



qucncy is reduced towards the zero mechanical excitation 
limit. One would expect the same to happen for the het- 
erogeneities of the displacements. However, we have also 
seen that the system size limits the largest heterogeneities 
of the displacements. And indeed, when we reduce the 
vibration, there is saturation of the displacement hetero- 
geneities, while those for the contacts increase fourfold 
(figure [IHi;c) and (d)). 
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(a) (b) 

FIG. 19: Dynamical correlations of the contacts at 
short and long times, (color online) (a): Dynamic sus- 
ceptibility of the contacts xl vs. the lag time r. Same pack- 
ing fractions as in figure [S] (b) : Dynamic susceptibility at 
short time xi" (x) and maximal dynamical susceptibility xl* 
(□) vs. reduced packing fraction, e. The vibration frequency 
/ = 10 Hz, i.e. 7 = 1.4. 



C. Short time origin of the heterogeneities 

Figure llTf a) indicates that the nonmonotonic depen- 
dence of X4 on the packing fraction applies not only at 
timescales of ^ r*" , but is also manifested at the short- 
est timescales of the data acquired stroboscopically, that 
is for ~ one cycle or 10^ microscopic times. The same 
holds true for the contacts. Figure [TOla) and (b) re- 
spectively display xH't) and xl^ = xli^o): together with 
xl* = xH't'^*) functions of the packing fraction: both 
are nonmonotonic, suggesting that the dynamical het- 
erogeneities of the contact dynamics have roots in the 
structure of the contact network. Still, xt^ is smaller 
than xl*, indicating that the heterogeneities, present at 
short time, build up progressively via a process which 
remains to be explained. 

The above results suggest that the contact network 
itself is heterogeneous. Whereas a number of papers dis- 
cuss the heterogeneities of the force network in terms of 
the force intensities, we are not aware of a detailed exami- 
nation of the spatial correlations in the contact network. 
A map of the instantaneous contact network for pack- 
ing fractions lower than cj)^ is provided in figure I20f a) . 
One immediately notices a rather heterogeneous organi- 
zation, with rather large holes where there are very few 
contacts. After interpolating the contact number on a 
grid, we compute the radial depcndance of its spatial au- 
tocorrelation G|(r). This quantity decays exponentially 
towards zero (figure I^DTb)). with a typical decay length 
^1 defined as 02(^1) = 0.2. ^| is non-monotonic with re- 
spect to packing fraction (figure l^c)). and has a smaU 
maximum at e* : the spatial correlations of the contacts 
are maximal at e*. An alternative and stronger evidence 
is provided by the static susceptibility 

N _ , , 

^2 = ~n — N v ^^^*^*' ^ ' 
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FIG. 20: Spatial correlations of the contacts, (color 
online) (a): Instantaneous map of the contact number, for 
e = —0.0091. The color map varies from white {zi{t) = 0) to 
black {zi{t) > 6). Contacts links are indicated in red. (b): 
Spatial correlations of the contacts Gi vs. r. Inset: Gi vs. 

Same packing fractions as in figure [S] (c): Spatial cor- 
relation length of the contacts vs. reduced packing fraction 
e. (d): Contact susceptibility X2 vs. reduced fraction e. The 
vibration frequency / = 10 Hz, i.e. 7 = 1.4. 



where var z{t) is the temporal variance of the instanta- 
neous average number of contacts, z{t) = -i- ^i(^)- 
X2 is again maximum at the packing fraction e* {&g- 
ure YIOi d]]. pointing at a maximum static correlation. 

Altogether, the present results: first, confirm that the 
dynamical heterogeneities observed in the displacement 
fields are connected to the heterogeneous dynamics of the 
contact; second, they indicate that the heterogeneities 
are already present in the static properties of the contact 
network. Such a connection is quite remarkable, and it 
would be interesting to see whether a similar one exists in 
the case of thermal soft spheres close to Jamming. Also, 
the mechanism by which the static behavior at short time 
and the dynamics at longer time are connected remains 
unclear, and would deserve further investigation. 



VII. DISCUSSION 

We recall here the motivations which lead us to con- 
duct this comprehensive study of vibrated photo-clastic 
disks. 

First we sought to confirm our first observations of 
dynamical heterogeneities in a very dense system of vi- 
brated brass disks [23l - [25j . These heterogeneities are 



rather singular in the sense that they concern very small 
displacements, of the order of 10~^ grain diameters, and 
occur for very large packing fraction as compared to those 
observed in other granular systems [s^ . There ex- 
ist no other experimental evidence of such dynamical 
heterogeneities, except perhaps in one colloidal experi- 
ment [50[, and in other quasi-static experiments by the 
authors elsewhere [5l| . but it remains unclear whether 
these different experiments probe the same physics. Only 
recently [l|, \^ , similar observations have been reported 
in numerical simulations of soft spheres, a system that is 
a priori rather different from vibrated granular media. 

Our first set of experiments conducted with soft photo- 
elastic disks confirmed results for a system consisting 
brass disks, and lead to the observation that the simi- 
larities with the simulations of thermal soft spheres were 
stronger than expected [l| . Since the authors of that nu- 
merical study argued that existing colloidal experiments 
are rather far from the critical regime of Jamming, ei- 
ther because the packing fractions are too loose, or be- 
cause the temperature is too high, we chose to decrease 
the vibration frequency in our system and to explore the 
vicinity of the zero excitation limit. Indeed, one would 
like to know to what extent thermal harmonic spheres 
have anything to say about the dynamical criticality of 
the granular packings, and conversely, whether granular 
experiments can provide physical insight into the ideal 
system of harmonic spheres. 

The discussion section below is organized as followed. 
After a brief synthesis of the results, we first compare and 
reconcile the observations performed for the hard (brass) 
and the soft (photo-elastic disks), before discussing the 
analogy between the thermal soft spheres and our exper- 
imental systems. 



A. Synthesis 

We have conducted systematic experiments of horizon- 
tally vibrated grains, decreasing the packing fraction over 
a very small range of high packing fractions where the dy- 
namics of both the contacts and of the positions of the 
grains is frozen. Despite a strongly anisotropic mechani- 
cal forcing at large scales, the system at the scale of the 
grain is isotropic: nonlinear mechanisms, together with 
disorder, redistribute the energy at small scales, causing 
the system to progressively lose any memory of the forc- 
ing anisotropy. This is roughly analogous to the energy 
cascade in turbulence. 

As previously noted, by using fast stroboscopic acqui- 
sition, we computed the average displacements over more 
than six temporal decades, once the short term oscillating 
dynamics and long term convection have been removed. 
We clearly identify a ballistic regime, followed by a long 
plateau, eventually followed by a crossover to a very long 
time diffusive regime for low enough packing fractions. 
These observations allowed us to measure the size of the 
cages. A, as a function of the packing fraction in several 
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independent ways. 

Within the timescales where the grains are trapped in 
their cage, two distinct crossovers are observed. One is 
"structural" in the sense that it is revealed by the aver- 
age number of contacts, which starts increasing sharply 
at the packing fraction cf)^ . The second is "dynamical" 
in the sense that it is indicated by a maximum of the 
dynamical heterogeneities of both the contacts and the 
displacements at a packing fraction (jf < 0^. We have 
demonstrated that the "dynamical" crossover is rooted 
first in the structure of the contact network, and sec- 
ond, that it is related to the spatial fluctuations of the 
contact number. By contrast, the "structural" crossover 
is given by its average value. Both signatures converge 
to a unique packing fraction when the excitation is re- 
duced towards the zero excitation limit. We interpret 
this packing fraction as the Jamming transition for the 
present experimental system and compression protocol. 
The critical nature of the transition is suggested by the 
sharp increase of the dynamical susceptibilities when the 
vibration is reduced towards the zero excitation limit. 
The two crossovers can be seen as the analogs of the 
Widom lines reported in the supercritical region of equi- 
librium phase transitions (26l |44 | . 



B. Soft vs. Hard 

In an earlier experiment, within the same apparatus 
but with hard (brass) disks [2^, [l^l, the authors re- 
ported the first experimental evidence of dynamical het- 
erogeneities involving very small displacements of grains, 
within a structure almost completely frozen. These dy- 
namical heterogeneities were rather different from those 
observed close to the glass transition, and the authors 
correctly attributed their observation to Jamming. How- 
ever, they could not precisely identify the underlying 



mechanism responsible for these heterogeneities. The 
present study has clearly demonstrated that the hetero- 
geneities have their origin in the dynamics of the contact 
network. Also, the existence of this maximum in dynam- 
ical heterogeneities suggests that the experiment probed 
both sides of the Jamming transition, a puzzling conclu- 
sion given the stiffness of the brass disks. The present 
study with soft disks solves this apparent contradiction 
in the following way. We have seen that there are several 
signatures of point J at finite mechanical excitation, 7, 
and that the one associated with the dynamical hetero- 
geneities occurs at a lower packing fraction, 0*(7), than 
the one at which the average number of contact increases, 
0t(7). 

In the case of the brass disks, the authors reported (see 
figure ETT a)) that the maximum of the dynamical hetero- 
geneities occurs for the packing fraction where Pdyn (0) 
and Pstat{4') intersect. This is also the case for the 
soft disks (see figure [HJb)): the experiment with the 
brass disks actually probed the dynamical crossover (f>*, 
both sides of which lie below the structural signature of 
the Jamming transition. In the case of brass disks, it 
is not possible to measure the average number of con- 
tacts. However, assuming Hertz' law, the stiffness of two 
compressed 4 mm height cylinders made of brass (Young 
modulus, E = 100 GPa) is hrass ~ 3 x 10^ N/m. By 
comparison, the stiffness of the force sensor and piston 
system is kpiston ^ 6.10^ N/m and the brass grains can 
be considered as hard. In that case. Jamming is the point 
at which the pressure diverges [ij, [13 , and the packing 
fraction at which the pressure sharply increases (sec iig- 
ure l^lTa)). provides a good estimate of the structural 
crossover 0^ . 

We also note that the range of packing fractions over 
which the crossovers are observed are very different. The 
crossovers occur for lower packing fractions and on a 
broader range in the case of the soft disks than in the 
case of the hard particles. This is not so surprising, given 
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(a) (b) 

FIG. 21: Hard vs. Soft, (color online) Piston force (top) 
and Maximal dynamical susceptibility of the displacements 
(bottom) vs. reduced packing fraction, e, for (a): hard brass 
disks [2J] and (b): soft photo-elastic disks. (O)- Ptot, (□): 
PsTAT, (A): Pdyn as in figure [G] The vibration frequency 
/ = 10 Hz, i.e. 7 = 1.4. Dashed lines indicate e* and e = 0. 
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(a) (b) 

FIG. 22: Towards zero vibration (color online) (a): 
Maximum of the dynamical susceptibility of the contact 
maximum(x4*) for soft grains (x) and hard grains (0)i esti- 
mated from max(x4 )/20, versus the split |e*| between static 
and dynamics signatures of Jamming, (b): MSD Plateau vs. 
density e, for 7 = 0.5 (Q), 7 = 0.8 (x), 7 = 1.4 (□) and for 
hard brass disks at 7 = 1.4 (O). 
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that the friction coefficient between the grains and be- 
tween the grains and the glass plate are different. The 
soft disks have indeed a larger friction coefficient, so that 
their Jamming transition is expected for lower values of 
the packing fraction [H, [H^. They also have a larger 
friction coefficient with the glass plate shaking them so 
that the energy transfer and dissipation are different. It 
is remarkable that, despite these differences, when plot- 
ting the peak of the maximal dynamical susceptibility 
of the displacements as a function of the split separat- 
ing the dynamical and the structural crossovers (see fig- 
ure im^a)), the experiment conducted with brass disks 
align with those conducted with the soft disks at differ- 
ent vibration frequencies. This suggests that the hard 
disks vibrated at a frequency / = 10 Hz, i.e. 7 = 1.4 be- 
have as soft disks with a much smaller effective value of 
7: the injection of energy is much less efficient in the case 
of the hard, less frictional, disks. It also indicates that 
friction plays a role in the absolute value of the packing 
fraction t/fj, as well as in the efficiency of the mechani- 
cal excitation, but not in the physics observed at finite 
vibration. 

We can further confirm that the hard disks behave like 
the soft disks at a lower level of excitation, by compar- 
ing the mean square displacement of the grains in 
the plateau regime, also called the Dcbye- Waller factor. 
For the brass disks, no fast camera acquisition were con- 
ducted, but one can take the displacements over one vi- 
bration cycle as an upper bound of the plateau value. 
Figure \2^h) displays this Debye- Waller factor for the 
three experiments using the soft photo-elastic disks at 
three different vibration frequency, as reported in the 
present paper, and for the experiment with the hard brass 
disks. The value of ^ 10~^ is significantly lower for 
the brass disks, confirming that they sit closer to the zero 
vibration limit. 



C. A-thermal vs. Thermal 

We have just seen that the physics of the Jamming 
transition of granular media in the presence of vibra- 
tion is robust with respect to the specific properties of 
the grains. However, as stated in the introduction, the 
Jamming transition is precisely defined and well charac- 
terized for thermal soft spheres, not for frictional grains. 
To what extent does it describe our present observations? 
In other words, docs the street lamp of thermal Jamming 
illuminate the granular world? 

To answer this question we return to the recent nu- 
merical simulations [l[ in which the authors studied the 
dynamical behavior of thermal soft spheres, close to Jam- 
ming, when approaching the limit of zero temperature. 
They report the existence of dynamical heterogeneities 
of the displacements when probed at very small scales. 
The peak of the maximal heterogeneities increases when 
the temperature is decreased to the T = limit, and 
the packing fraction at which this peak occurs decreases 



in the same limit. All these observations are identical to 
those reported in the present work. Unfortunately, in the 
experiments, we don't have access to a well defined and 
unique value of the Jamming packing fraction at strictly 
zero vibration: it varies from one realization to another. 
Also, comparing the scalings with the distance to point 
J would require a proper definition of an effective tem- 
perature, a notoriously difficult, if not impossible, task. 

However, we can follow Ikeda ct al. [ij, who use the 
mean square displacement in the plateau regime as a 
sensitive thermometer close to Jamming, and when they 
compare their observations to experimental colloidal sys- 
tems [ssl - ISTj . From A^ (0) computed for several tempera- 
ture, and knowing the range of packing fraction explored 
by the colloidal experiments, it is straightforward, using 
the mean square displacements reported for such exper- 
iments, to locate them in the Temperature-packing frac- 
tion parameter space. The advantage of such a method 
is that is does not require any information about the 
details of the interaction potential, nor the knowledge 
of the kinetic energy. We follow exactly the same pro- 
cedure. The authors report that the mean square dis- 
placement in the vicinity of Jamming decreases from 
10"'^ to 10^^ particle diameters when the temperature 
is decreased from T = IQ-^ to T = 10"®. As ob- 
served in figure [22l bl. for the soft photo-elastic disks 
experiments, A^ € [10~*10~^], corresponding to temper- 
atures of [10~^10~^] and, for the hard brass disks one, 
A^ ^ 10~^, corresponding to a temperature T 10~*. 
Of course, these temperatures have no thermodynamic 
meaning, they are essentially a measure of the kinetic 
energy at short times. Figure [13] summarizes the above 
discussion: the granular experiment indeed probes the 
critical regime of Jamming at finite temperature. 

One must realize the impact of such a conclusion. 
Shaken granular systems and thermal soft spheres are 
very different. In large part, due to the effect of dissipa- 
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FIG. 23: Temperature-density phase space A demon- 
stration that the present granular experiments do probe the 
criticality of the Jamming transition. Adapted from [l| . 
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tiong/friction, shaken granular media are out of equilib- 
rium systems for which detailed balance does not hold. 
Energy is injected at large scale and dissipated at small 
scales. In the present case, this ensures the isotropy of 
the displacements at short time, but it is also responsible 
for the large convection pattern that we have removed. 
Such effects would never have existed in an equilibrium 
glass. Despite these significant differences, the dynam- 
ics seems to obey the same physics as soon as a small 
amount of vibrations or thermal agitation is present. 

VIII. CONCLUSION 

In this work, we have demonstrated that, in the pres- 
ence of agitation, the Jamming transition's singular fea- 
tures are blurred into two crossovers, a structural one 
indicated by the increase of the contact number, which 
is directly inherited from the zero excitation limit case, 
and a dynamical one, specific to the presence of agitation. 
The contact network develops heterogeneous dynamics, 
which in turn induce heterogeneous displacements at very 
small scale. These heterogeneities take place within the 
vibrational regimes, while the structure of the glass re- 
mains essentially frozen, and they are related to struc- 
tural heterogeneities in the contacts network itself. 

These observations match the recent results reported in 
numerical simulations of harmonic soft spheres [l| very 



well, and we demonstrated that the critical regime of 
point J is indeed probed by our granular experiments. 
The street lamp of Jamming illuminates the granular 
world. This strongly suggests that similar experiments 
be conducted in other systems, such as foams or emul- 
sions, provided one finds a way to "vibrate" them. Note 
that in all cases, it is key to reach a spatial resolution of 
the order of a thousandth of the size of the elementary 
component. At present this limitation has for instance 
prevented colloidal experiment from probing the critical 
regime, although in principal, colloidal suspension are the 
closest systems to the thermal soft spheres. 

Finally, one may wonder whether similar conclusions 
apply in the presence of an external driving force. Con- 
ducting shear experiments in the limit of very weak vibra- 
tions, in the spirit of [ssj . while monitoring the contact 
dynamics is a promising perspective. 
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